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Abstract 

The intrinsic decay rate of orthopositronium formed in Si02 powder is measured 
using the direct 27 correction method such that the time dependence of the pick-off 
annihilation rate is precisely determined using high energy-resolution germanium 
detectors. As a systematic test, two different types of Si02 powder are used with 
consistent findings. The intrinsic decay rate of orthopositronium is found to be 
7.0396 ± 0.0012(stat.) ± 0.0011(sys.)/is~\ which is consistent with previous mea- 
surements using Si02 powder with about twice the accuracy. Results agree well with 
a recent 0{a'^) QED prediction, varying 3.8 — 5.6 experimental standard deviations 
from other measurements. 
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1 Introduction 



Positronium (Ps), the bound state of an electron and positron, is a pure 
quantum electro-dynamical system providing a highly sensitive field for test- 
ing bound state quantum electrodynamics (QED). The triplet (I'^Si) state 
of Ps, orthopositronium (o-Ps), will in most cases decay into three photons 
due to odd-parity under the C-transformation. Because o-Ps has a lifetime 
about 1140 times longer compared to the singlet state, parapositronium (p- 
Ps), this enables direct and precise measurement of the intrinsic decay rate of 
orthopositronium, Ao-ps [1,2], although obtained values are much larger, i.e., 
5.2 and 8.2 experimental standard deviations, than a recent non-relativistic 
QED calculation (7.039 979(11) iis'^) corrected up to 0{a^) [3]. To eluci- 
date discrepancies, a variety of experiments have since been carried out to 
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search for the exotic decay mode of o-Ps, resulting in no evidence so far 
[4,5,6,7,8,9,10,11,12,13]. 



As some fraction of o-Ps inevitably results in 'pick-off' annihilations due to 
collisions with atomic electrons of the target material, the observed o-Ps de- 
cay rate Xobs is a sum of the intrinsic o-Ps decay rate Aq-ps and the pick-off 
annihilation rate into 27's, Xpick, i-e., 

^obs{t) — Aa-y -|- Xpick{t). (1) 



Contributions from exotic decays are assumed to be zero because they are con- 
firmed to be less than 200 ppm [4,5,6,7,8,9,10,11,12,13]. Xpick{t) is proportional 
to the rate of o-Ps collisions with the target materials, i.e.; Xpick = iT'(^av{t), 
where n is the density of the target, aa is the annihilation cross section, and 
v{t) the time dependent velocity of o-Ps. Due to the thermalization process of 
o-Ps, this necessitates expressing Xpick as a function of time whose properties 
are dependent on the surrounding materials. 

In previous measurements [1,2], Aofes's were measured by varying the densities 
of the target materials (n) or size of the cavities. For the uniform distribution 
of o-Ps in the cavity, Xpick is depending on the size and shape of the cavity due 
to the relation, v{t) — {v{t)/A)S/V, where v{t) is the average o-Ps velocity, S 
is the cavity surface area, and V is the volume of the cavity. 

The extrapolation to zero gas density or infinite cavity size was expected 
to yield the decay rate in a vacuum, Aa-y, under the assumption of quick 
thermalization with constant o-Ps velocity. However, this assumption con- 
tains a serious systematic error as pointed out in reference [14] in that it 
would make the obtained decay rate larger than the true value; i.e., the 
time-spectrum fitting was carried out using a simple exponential function, 
dN{t)/dt — A^oexp(— As-y t) + C, where C represents the fiat spectrum of acci- 
dental events. In order to cleanly eliminate the contribution from the prompt 
annihilations and fast p-Ps decay component, the decay spectrum needs to be 
fitted for times after these components decay away. A somewhat fast, begin- 
ning time of the fitting range was adopted to obtain the high statistics: the 
cavity experiment typically used t start — 170 ns for all runs [2], while the gas 
experiment [1] used tgtart = 180 ns for all gases except the lowest pressure A^e 
and neopcntane runs. Because it takes several hundred nsec for o-Ps to ther- 
malize via elastic collisions with surrounding molecules to an ultimate energy 
of about 0.03 eV depending on the type of target material [15,16], such fast 
t start tends to result in larger values of X^^. Moreover, in the cavity experiment, 
o-Ps is not well thermalized due to the small escape and collision rates. 

The new technique called the direct 27 correction method was introduced 
to overcome the problem regarding the extrapolation. The following section 



2 



briefly explains its concept and procedure. Once a precise thermalization func- 
tion is obtained, \pick{t) wiU contain all information about the process. The 
population of o-Ps at time t, N{t) can be expressed as 



Ar(t) =7V^exp -A37 



/ 








(2) 



The energy distribution of photons from the 3-body decay is continuous be- 
low the steep edge at 511 keV, whereas the pick-off annihilation is 2-body 
which produces a 511 keV monochromatic peak. Energy and timing infor- 
mation are simultaneously measured with high-energy resolution germanium 
detectors such that Xpick{t) / X^y can be determined from the energy spectrum 
of the emitted photon. Providing the ratio is determined as a function of time, 
the intrinsic decay rate of o-Ps, Xs^, can be directly obtained by fitting the 
observed time spectrum. The direct 27 correction method makes any extrap- 
olation procedure unnecessary and precluding use of the hypothesis of linear 
dependence at vacuum limit or o-Ps behavior at the material surface. 



2 Experiment 

The direct 27 correction method first applied in 1995 [17] provided consistent 
results with QED predictions, yet several problems remained: (i) accuracy 
was 400 ppm, being worse than those of the other experiments [1,2], (ii) de- 
cay rates systematicafly increased before tf^tart = 200 ns due to an unknown 
reason, and therefore to remove this uncertainty, final results were obtained 
using data after 220 ns, and (ii) systematic error regarding the Stark effect 
was not estimated. Improving the method by considering these problems will 
be described in the later sections. We employed a ^^Na positron source (dia., 
2 mm) with approximate strength of 0.03;LiCi, being sandwiched between two 
sheets of plastic scintillators (NE104) and held by a cone made of aluminized 
mylar (Fig. 1). The scintillators and mylar were 100 /im- and 25 /^m-thick, 
respectively, corresponding to a total area density of about 17 mg/cm^. The 
cone is situated at the center of a cylindrical 50 mm-diameter vacuum con- 
tainer made of 1 mm-thick glass, being filled with Si02 powder and evacuated 
down to 5 X 10^^ Torr. Two different types of Si02 powder (Nippon Aerosil 
Ltd.) were used (Table 1), with the biggest difference being in the mean dis- 
tance between grains (340 vs. 200 nm) such that different pick-off ratio would 
be obtained. Using these powders, two 6-month runs were performed (RUNl 
and 2). To remove water molecules absorbed on the grain surface, the powder 
was heated just before evacuation. 

While travehng through the scintillators, the source-emitted positrons deposit 
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p =0.03 g/cm^ for RUN I 
0.05 g/cm^for RUN I 



Fig. 1. Schematic diagram of experimental setup. D, t represent the diameter and 
the thickness (longitudinal length) of each cylindrical component respectively. The 
density, p, and the other characteristics of the silica powder are summarized in 
Table 1. Upper right figure illustrates the geometrical location of each component 
viewed from top. 



RUNl RUN2 

primary grain size (nm) 
surface area {m?/g) 
density (g/cm^) 
mean distance between grains (nm) 
surface 

Table 1 

Characteristics of Si02 powders used in the measurements. 

energy from 50 — 100 keV which produces scintillation photons transmitted 
to the trigger PMT (Hamamatsu H3165-04). The magnetic field around the 
assembly is measured to be 0.5±0.1 Gauss, equivalent to terrestrial magnetism 
which barely contributes to the mixing between o-Ps and p-Ps estimated to 
be 3 X 10~^^. The observed o-Ps is therefore regarded as a pure o-Ps sample. 



7 7 

300 ± 30 260 ± 30 

0.03 0.05 

340 200 

hydrophile hydrophobe 



Two high-purity coaxial germanium detectors (Ortec GEM38195), referred to 
as GeO and Gel, precisely measured the thermalization process, respectively 
having crystal sizes of 58.3 and 60.3 mm in dia. and 73.8 and 67.4 mm in length, 
and a solid angle of 5.7 and 6.0% of An. Energy resolutions were measured 
using several line 7 sources, with typical resultant values of 0.53 and 0.64 keV 
in sigma for GeO and Gel, respectively, at 514 kcV. Lead sheets 2.0-mm- 
thick were placed in front of each detector for suppressing contributions from 
simultaneous low-energy 7 hits from the 37 decay of o-Ps. 

Four large cylindrical Nal(Tl) scintillators (Scionix 127A127/5; 127 mm (0) x 
127 mm {t)) simultaneously measured the time and energy information from 
each decay. Due to their higher efficiency and faster time response, the time 
spectrum can be fit to the o-Ps decay rate determination. At 662 keV, their 
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Fig. 2. Simplified sclicmc of the trigger system. In figure, DIS and LOGIC represent 
the discriminator and the logic unit components respectively. Descriptions of signal 
flows are found in the text. 



energy resolution is typically 22.0, 24.4, 21.4, and 32.0 keV. 

A new time- to- digital converter (TDC) was employed, jointly developed with 
the High Energy Accelerator Research Organization (KEK) and customized 
for the present measurement. This direct clock-counting type TDC, referred to 
as the 2-GHz TDC, has a time-resolution of 0.5 nsec with known accuracy of 
l-ppm, equipped with 8 channels being able to accommodate multiple signals. 
The time range for each channel is 32 fis, and the integral non-linearity (INL) 
is expected to be extremely small at < 15 ppm. To provide a systematic check, 
we used a 200-MHz internal clock-based TDC (Hoshin C006) with 5 ns time 
resolution. 



The electronic system consists of three sections: the trigger, Nal(Tl), and 
Ge sections. Each has a latch veto structure to prevent overlapping events. 
The main latch of the trigger section is set by the signals themselves and 
released with system reset. To ensure precise measurement, the latches of the 
Nal(Tl) and Ge sections are also set by themselves and released by reset signals 
from the trigger section. The trigger PMT signal is fed into a fast leading 
discriminator whose output provides common start signals for the TDCs as 
well as their reset signals. Each output signal from the photo-multipliers to 
the NaI(Tl)s is fed into three ADCs and a fast leading edge discriminator 
that provides stop signals for the TDCs and gate timing for the ADCs. One 
ADC with a S-fxs gate width, called the NAI-WD (Lecroy 2249W), measures 
the whole charge for the duration of the signal, with slewing correction being 
carried out using energy information. The other two ADCs eliminate pile-up 
events at the tail of the signal and base-line fluctuations, respectively. The 
NAI-NW (mainly Lecroy 2249W) ADC measures the signal charge with a 
narrower gate width (250 ns for RUNl and 400 ns for RUN2) comparable to 
the intrinsic decay time constant of the Nal(Tl) scintillator, while the NAI-BS 
(Lecroy 2249A) ADC measures the base-line condition of the signal (180 ns 
width) just prior to the event. All gates for these ADCs are individual and 
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gate timings are optimized for signal timings. 

Signal outputs from the Ge detectors are used for timing measurements and 
obtaining precise energy information. One signal is fed into a fast-filter ampli- 
fier (FFA)(Ortcc 579) whose output is used as the stop signal for the TDCs, 
two auxiliary ADC signals, and three discriminators whose thresholds are 
set at —50, —100, and —150 mV. Each FFA output is fed into another high- 
resolution TDC having fine resolution of 250 ps. Signal timings are determined 
utilizing signal shapes calculated from different threshold discriminating times 
and extrapolating to the intrinsic timing [18] . Good time resolution of 4 ns is 
obtained, as is efficient rejection power for slow rise signal components known 
to disturb Gc timing spectroscopy. Similar to the Nal(Tl) detectors, the two 
auxiliary ADCs effectively reject pile-up events by measuring the signal charge 
with a narrower gate (2.0(1.3) jjs width for RUN2(1)) and earher timing gate 
(200 ns width just prior to the signal). Precise time and energy resolutions 
can be obtained with this rejection scheme. The other output from each Ge 
detector is independently fed into a SHAPER spectroscopy amplifier (Ortcc 
673). Amplified with a 6-iis time constant, the uni-polar output is provided 
to a Wilkinson type peak-holding ADC called GE-WD that provides a precise 
energy spectrum. For RUNl and RUN2, 2.7 x 10^ p+ and 2.0 x lO'^ /3+ events 
were recorded, respectively Each run had a total data acquisition period of 
about six months. Area temperature was maintained within ±0.5°C to ensure 
stability of the amplifiers, ADCs, and TDCs. 



3 Analysis and discussion 

The ratio \pick{t) / ^S'y is determined using the energy spectrum measured by 
the Ge detectors. The energy spectrum of the o-Ps decay sample, referred to 
as the o-Ps spectrum, is obtained by subtracting accidental contributions from 
the measured spectrum. The 37-decay continuum spectrum is calculated using 
Monte Carlo simulation in which the geometry and various material distribu- 
tions are reproduced in detail. Three photons are generated according to the 
0{a) corrected energy spectrum [19], followed by successive photo-electric, 
incoherent, and coherent interactions with materials until all photon energy is 
either deposited or escapes from the detectors. The response function of the 
detectors is determined based on the measured spectrum of monochromatic 
7-rays emitted from ^^^Eu, ^^Sr, and ^^^Cs, with this function being used in 
the simulation. We refer to the obtained spectrum as the S'y-spectrum which 
is normalized to the o-Ps spectrum with the ratio of event numbers within the 
energy domain 480-505 keV. The region is chosen to minimize, the Compton 
scattering effect from pick-off events, dependence of the detector efficiencies 
on energy, and the 0{a) correction to the 37-spectrum. Figure 3(a) shows 
good agreement between the o-Ps spectrum and 37 spectrum below 508 keV, 
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Fig. 3. (a) Energy spectrum of o-Ps decay 7's obtained by Ge detectors for RUNl. 
Dots with error bars represent data points in a time window of 150—700 ns. The sohd 
line shows the 37-decay spectrum calculated by Monte Carlo simulation. Shaded 
area indicate simultaneous hits of two low-energy 7's from a 8-7 decay estimated 
using the same simulation, (b) Pick-off spectrum obtained after subtracting the 87 
contribution from the o-Ps spectrum for RUNl. The solid line represents the fit 
result. 



where the pick-off annihilation peak is evident at the edge of the 37-decay 
continuum. It should be noted that the lead sheets in front of the Ge detec- 
tors effectively suppress the contribution from events in which two low-energy 
7's emitted from a 37-decay simultaneously hit one detector and that the 
37-spectrum well reproduces such simultaneous events. 

Figure 3(b) shows an enlarged view of the observed o-Ps spectrum after sub- 
tracting the 37-spectrum, where good agreement with a detector response 
function [20,21] is present. The centroids of the peak spectrum are obtained 
as 510.997 t^^^ (keV) for RUNl and 510.995 1^^^ (keV) for RUN2. Good 
consistency with 511.0 (keV) indicates successful subtraction of the 37 con- 
tribution such that the resultant peak can be regarded as pure pick-off an- 
nihilation samples. Obtained ratios of Xpick/^s-y at a typical time window of 
150-700 ns are 0.01049(8) and 0.01582(9) for RUNl and RUN2, respectively. 
The weighted mean time within the range for each run is 287.5 ns (RUNl) 
and 284.3 ns (RUN2). The large difference between the two runs originates 
from a different mean distance of powder grains (Table 1) since the lambda 
ratio is a function of time. 

Figure 4 shows results when the procedures described above for a typical time 
window were applied to multiple, narrower time windows. Since the fractional 
energy loss of o-Ps per collision with Si02 powder and the collision rate are 
both dependent on its energy, the time dependence of the average kinetic 
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energy of o-Ps at time t, E{t) can be derived from the Boltzmann equation 



^^Eit) = -^2mpMt) (^E{t) - ^knT^ £ a, (^^J^ , (3) 

where mps is the mass of o-Ps, T room temperature, and kp the Bohzmann 
constant. The last term, the momentum transfer cross section of Si02 is ex- 
panded in terms of E{t), i.e., the coefficients aj represent the effect of effective 
mass at the surface of the Si02 grain and mean distance between the grains. 
Since the pick-off rate is proportional to the average velocity of the o-Ps, the 
ratio 9{t) = Xpickit) / ^s-y can be expressed by a differential equation [22,23], 
i.e., 

j^e{t) = -c (9{tf - 91) e{tr, (4) 

where C is a constant, 6oo = 0{t ^ oo), and the last summation term in Eq. 3 
is replaced with an arbitrary real number /?. The measured \pickit) /Xs-y^s are 
fit with Eq. 4, and the first-order differential equation is numerically solved 
using the Runge-Kutta method. The data point at a time window of 40 — 45 ns 
is used for a fixed initial value. Fitting is carried out using four parameters: 
C, 9ao, (3, and the deviation from fixed point 6 found to be —0.27 x 10"^ and 
—0.32 X 10^"^ for RUNl and RUN2, respectively. These small values indicate 
that the fitted values rarely depend on the initial condition of the Runge-Kutta 
method. 

Figure 4 shows best fit results using the MINUIT library [24], where the pick- 
off rate cannot be assumed as constant, even in Si02 powder where the collision 
rate is expected to be higher. Table 2 gives resultant values for these fitting 
parameters. 



Parameters 


RUN 1 


RUN 2 


C 

2/3 
xVd.o.f 


(0.166 X 10-2 
(0.761 1°:°?°) X 10-2 
2 64 +°-^^ 

^■^^ -0.30 

9.53/(23-4) 


(0.522 10-[|}6) X 10-2 
(0.999 l°i2) X 10-2 

9 -I 9 +0.25 
^■^■^ -0.22 

12.27/(23-4) 



Table 2 

Summary of Xpick{t) / X^-y fitting results. 

Figure 5 shows the time spectrum of Nal(Tl) scintillators for RUN 1 with 
an energy window of 370 — 440 keV, where a sharp peak in prompt positron 
annihilation is followed by the exponential decay curve of o-Ps and then the 
constant accidental spectrum. The o-Ps curve is widely observed over about 
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time (ns) 



Fig. 4. The ratio Xpick{t) / X^^ are plotted as a function of time. Open circles arc data 
points for RUNl and closed circles for RUN2. Solid lines represent best fit results 
obtained using Eq. 4. 
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time (ns) 

Fig. 5. Time spectrum of Nal(Tl) scintillators for RUN 1 within an energy window 
of 370 — 440 keV. The inset shows an enlarged view of the prompt peak, where good 
time resolution of o" = 2.2 ns is obtained. 

1.2 us corresponding to about eight times the o-Ps lifetime, being partly due 
to the use of a weak positron source (0.03 nCi) and good suppression of 
accidental contributions by selecting the 7-energy. To effectively eliminate 
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Fig. 6. Decay rates as a function of fitting start time. Small horizontal lines on the 
error bars represent the size of statistical errors solely due to fitting, while vertically 
extending bars include propagated errors from Xpick{t) / determination. Open 
and closed data points indicate values for RUNl and RUN2, respectively. Dashed 
and solid lines show one standard deviation obtained at t = 100 ns for RUNl and 
RUN2, respectively. 

pile-up events, a base-line cut condition was applied, and to further reject 
pile-up events, events with small differences between two ADC values (wide 
and narrow gates) were selected. 

After selecting events in the energy window of 370 — 440 keV, we fit resultant 
time spectrum using the least square method, i.e.. 



Nobs{t) = exp{-Rstopt) 



^37 



A 



37 



(5) 



where epick and e^^ are the detection efficiencies for pick-off annihilations and 
37 decays. Rgtop in the function is the measured stop rate of the Nal(Tl) 
scintillators, representing the fact that time interval measurement accepts only 
the first 7 as a stop signal. Xpick/^sj is about 1% due to the low-density of the 
Si02 powder, i.e., the ratio of error propagation to decay rate is suppressed by a 
factor of 100. The value of epick/ ^^-y is estimated using Monte Carlo simulation, 
yielding 0.0651(7) and 0.0607(6) for RUN 1 and RUN2, respectively. The time 
spectrum is corrected in terms of the Rstop (order of 1 kHz) before fitting. 
Fitting stop time is fixed at 5000 ns. 

Figure 6 shows obtained fitting as a function of fitting start time for both 
runs, where values are stable with respect to fitting start time except those 
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before 80 ns. Since the fitting x^'s for both runs rapidly increase before 90 ns 
due to the tail effect of prompt events, values at t = 100 ns are taken as the 
final results. The reduced x^'s at t = 100 ns are 1.006 and 1.051 for RUN 1 
and RUN 2, respectively. Note the good correspondence between runs. Since 
no large difference in thermalization properties are apparent in Fig. 6, this 
indicates that our method correctly takes into account thermalization. The 
obtained decay rates are As^ = 7.03991 ± 0.0017{stat.) iJ,s~^ for RUN 1 and 
As^ = 7.03935 ± Q.QQ17{stat.) /xs"^ for RUN 2. 



4 Discussion of systematic errors 



oOUrCc OI v^OIlTjilDUTjlOIla 


RTTl\n (•nryrr)\ 


riUlNZ yppnij 


TDC module dependence 






- Calibration 


< 1 


< 1 


- Stability 


2 ~ 3 


2 ~ 3 


— Integral Non Linearity 


< 15 


< 15 


- Differential Non Linearity 


Negligible 


Negligible 


Cut condition dependence 






- Base Line Selection 


-17 and +89 


-11 and +23 


- WD-NW condition 


-6 and +45 


-20 


Monte Carlo dependence 






— Normalization 


±99 


±113 


- Relative efficiency 






of Nal(Tl) scintillator 


+7 


+7 


- Inhomogeneity of Si02 powder 


< ±55 


< ±88 


Other Sources 






— Zeeman effect 


-5 


-5 


- Three-photon annihilation 


-26 


-42 


- Stark effect 


-21 


-6 


- n = 2 excited state 


+19 


+19 


Total 


-120 and +153 


-151 and +148 



Table 3 

Summary of systematic errors. See text for details of each item. 

The sources of systematic errors are summarized in Table 3 and each item is 
depicted in the following discussions. 

TDC module related errors: Calibration accuracy (< 1 ppm) and stability of 
the calibration oscillator (2 ~ 3 ppm) are known as they are product specifica- 
tions, while integral (< 15 ppm) and differential (negligible) non linearities are 
evaluated using basic performance of the 2-GHz TDC. Cut condition errors: 
Systematic errors for base-line selections and pile- up rejections were assessed 
by changing the cut criteria. The base-line cut was intentionally changed from 
±3(7 to ±5cr where the a represents standard deviation of the base-line dis- 
tribution. The decay rate results are compared with the nominal result with 
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±4(7, and systematic error is estimated as —17 ppm, +89 ppm for RUNl and 
—11 ppm, +23 ppm for RUN2. Regarding the pile-up events rejection, the 
WD-NW cut condition was changed from ±25keV to ±40keV to estimate de- 
viation from the nominal condition which is ±30 keV, and is estimated as 
—6 ppm, +45 ppm for RUNl and —20 ppm for RUN2. The dependence on 
the energy cut condition and trigger PMT threshold are considered to have a 
negligible effect. 

Monte Carlo simulation related errors: The predominant contribution to total 
systematic error is produced by uncertain normalization. That is, the number 
of pick-off events are determined by subtracting the normalized 8-7 spectrum 
of Monte Carlo simulation from the o-Ps spectrum, where changing the nor- 
malization factor affects the \ickif)/^z-i values and eventually propagates to 
the final result. Since the sharp fall-off of the 37-spectrum at 511 keV is solely 
produced by the good Ge energy resolution of cr = 0.5 keV, this subtraction 
only affects the lower side of the pick-off spectrum such that improper sub- 
traction results in asymmetry of the pick-off spectrum shape. Comparison of 
the asymmetries of the pick-off peak shape and the prompt peak annihilation 
spectrum is a good parameter for estimating this systematic error. The 1 cr 
error is assessed as ±99 ppm for RUNl and ±113 ppm for RUN2. 

Other MC simulation related errors: The relative efficiency of the Nal(Tl) 
scintillators. Spick/ has an uncertainty of +5% evaluated using a comparison 
of the energy spectrum of data and that of the simulation which reproduces 
data within an uncertainty of 2% except for the Compton free region where 
the deviation is 5% (+7 ppm). Si02 powder density in the MC simulation 
is conservatively changed by ±10% although the uniformity is known to be 
within a few %; a change resulting in an error of < ±55 ppm. 

Stark shift errors: The Stark shift stretches the hfetime of Ps atoms, i.e., a 
perturbative calculation shows that the shift is proportional to a square of 
the effective electric field E such that AAs^/Ag^ = 248 ■ {E/Eo f [25], where 
Eq = m^e^/Ii^ ~ 5.14 X 10^ V/cm. E is defined as the root-mean-square 
electric field sensed by o-Ps during its lifetime. Calculations have estimated 
two contributions exist based on measurements of the electrical charge-up 
on the primary grains of silica powders and electrical dipole moment on the 
surface of grains [18]. The charge-up is partly intrinsic depending on powder 
specifications and partly due to positron depositions from the (3'^ source. The 
effect, however, is negligible in both cases, i.e., on the level of 10~^ ppm at 
most. Silanol functional groups on the surface of the powder grain behave 
as an electrical dipole moment creating an effective field around the grains. 
Average densities are known to be as 2.5/nm^ and 0A4/nm^ for RUN 1 and 2, 
respectively. Accordingly, the effective field can be analytically calculated such 
that the contribution to the o-Ps decay rate is determined to be —21 ppm for 
RUNl and —6 ppm for RUN2. These estimations were confirmed using results 
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Fig. 7. Historical plot of o-Ps decay rate measurements including present results. 
Vertical line shows the 0(a^ )-corrected NRQED prediction [3]. Small vertical lines 
on the error bars indicate the size of errors solely due to statistics while error bars 
represent total ambiguities including systematic errors. 

from precise hyper-fine-structure (HFS) interval measurements of ground state 
Ps in silica powder [26], where the interval is proportional to the size of Stark 
effect. Considering the difference in powder densities used, the HFS results 
are consistent with our estimation. 



Other sources of systematic errors: Error contribution due to the Zeeman effect 
is estimated using the measured absolute magnetic field around the positro- 
nium assembly (—5 ppm). Since the 8-7 pick-off process can only occur at a 
certain ratio, the calculated relative frequency o"3^/o"27 ~ 1/378 is consistent 
with previous measurements [27], being —26 ppm for RUNl and —42 ppm for 
RUN2. The probability of the excited state (n=2) of Ps is about 3 x 10"^ [28], 
which could make the intrinsic decay rate 19 ppm smaller due to the low decay 
rate (|A3^) 

The above discussed systematic errors are regarded as independent contribu- 
tions such that the total systematic error can be calculated as their quadratic 
sum, resulting in —120 ppm, +153 ppm for RUNl and —151 ppm, +148 ppm 
for RUN2. 



5 Conclusions 



The decay rate of o-Ps formed in Si02 was measured using a direct 27 correc- 
tion method in which the thermalization effect of o-Ps is accounted for and in- 
tegrated into the time spectrum fitting procedure. Results were obtained using 
two runs with different types of Si02 powders, i.e., Ao-Ps(RUNl) — 7.03991 ± 
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0.0017(stat.)io:ooo8(sys.) fis-^ and Ao-ps(RUN2) = 7.03935 ± 0.0017{stat.) 
-o!ooii{sys.) jjis^^. Based on the results of both runs, a weighted average gave 
Ao-Ps = 7.0396 ±0.0012(stat.) ±0.0011(s7/s.)/is~^, being 1.8 times more accu- 
rate than previous Si02 measurement [17] . While this value agrees well with 
previous result, it disagrees with recent high-precision measurements, i.e., by 
5.6 a [1] and 3.8 a [2]. As illustrated in Fig. 7, our value agrees well with 
the NRQED prediction corrected up to O(a^) term [3]. In fact, the resul- 
tant decay rates were confirmed to be independent of fitting start time after 
t start — 90(80) ns for RUN1(2), and systematic increase in decay rate observed 
in previous measurement before tgtart — 200 ns [17] was eliminated using the 
new system; this improvement directly contributed to reduce the statistical 
error. 

Sincere gratitude is extended to Professor Toshio Hyodo, Dr. Yasuyuki Na- 
gashima, and Dr. Haruo Saito for their valuable suggestions and information 

regarding o-Ps properties and thermalization processes. Special appreciation 
goes to Dr. Masami Chiba for beneficial discussions on the o-Ps decay prob- 
lem and to Drs. Masahiro Ikeno and Osamu Sasaki for significant measurement 
improvement using the 2-GHz TDC module. 

Appendix 

During the submission procedure of this paper, we became aware of the recent 
result using nanoporous silica film by Ann Arbor group [29]. They obtained 
the value, Ao-ps = 7.0404 ± 0.0010(stat.) ± 0.0008(sys.) ns'^, which shows 
complete agreement with our result. 
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